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We investigate the expected neutrino emissivity from nuclear collisions in a magnetically domi- 
nated model of gamma-ray bursts motivated by recent observational and theoretical developments. 
The results indicate that small multi-GeV neutrino fluxes are expected for model parameter values 
which are typical of electromagnetically detected bursts. We show that for detecting at least one 
muon event in Icecube and its Deep Core sub-array, a single burst must be near the high end of 
the luminosity function, and at low redshifts z < 0.1, where the burst rate is very low. We also 
calculate the luminosity and distance ranges that can generate 0.01 — 1 muon events per GRB in 
the same detectors, which may be of interest if simultaneously detected electromagnetically, or if 
measured with future extensions of Icecube or other neutrino detectors with larger effective volume 
and better sensitivity. 



I. INTRODUCTION 

Observations of multi-GeV photons from gamma-ray 
bursts (GRB) recently accumulated by the Fermi satel- 
lite (e.g. [HIS]) have pointed out the need to re-evaluate 
the type of models used to explain the prompt pho- 
ton emission mechanisms and the location of the emis- 
sion regions in these objects (e.g. |3]). In particular, 
concerns about the radiative efficiency of usual internal 
shock models, and the larger radii required to avoid two- 
photon degradation of the spectra have spurred the inves- 
tigation of baryonic (non-MHD) jet models where the ra- 
diation arises in a jet photosphere In such baryon- 
loaded jet models, at a certain radius the timescale of 
the nuclear elastic collisions (which couple the proton p 
and neutron n components) becomes longer than the ex- 
pansion timescale, i.e. the collision optical depth falls 
below o(l), and the two components decouple from each 
other [ZllH]. The protons can continue to be accelerated 
by the radiation, while the neutrons, which have zero 
electric charge, start to coast with a constant Lorentz 
factor. Starting at this decoupling radius and for some 
distance beyond, the longitudinal drift velocity between 
the n and p becomes Aw > 0.5c, and they start to col- 
lide inelastically. Such large relative velocities between 
n and p components also can arise in realistic jets where 
the bulk Lorentz factor F depends on the polar angle 9, 
which also leads to inelastic collisions jHlini, as neutrons 
from the outer parts (sheath) of the jet thermally drift 
into the jet core. In both pictures, pions are created, 
which in the case of the dynamics being dominated by 
the baryons results in multi-GeV photons and neutrino 
production [51 |S]. 

A different approach towards resolving the radiative ef- 
ficiency of GRBs involves consideration of magnetically 
dominated jets |10H12| . Some of these magnetic models 
assume an almost baryon-free outflow jl2H16| . while in 



other cases a substantial but dynamically sub-dominant 
baryon load is assumed |101 [TT| fT7H22] . The baryons in 
such jet models are expected to accelerate at a differ- 
ent rate than in baryonic (non-MHD) jet models, and 
the different dynamics leads to quantitatively different 
predictions for the photon spectrum |23| . Similarly, it 
should lead to quantitatively different neutrino spectra, 
which we investigate in this paper. 

Unlike in the previous investigations cited above, here 
we consider the neutrino spectra arising from nuclear col- 
lision effects in magnetically dominated GRB outffows. 
In this case, both the radial np decoupling radius as 
well as the photosphere occur at larger radii from the 
central engine than in the non-magnetic case, and also 
the transverse drift of neutrons from the periphery of 
the jet into the jet core becomes important at different 
radii, where the physical conditions differ from those pre- 
viously considered. As a consequence, multi-GeV pho- 
tons are produced at somewhat softer energies and with 
appreciable time delays [53] respect to the MeV photo- 
spheric photons, but the detailed neutrino spectrum for 
such magnetically dominated jets has not been consid- 
ered so far. In this paper we investigate numerically the 
neutrino spectrum expected in magnetically dominated 
baryon-coUisional GRB models, taking into account both 
longitudinal and transverse n,p decoupling and inelastic 
collisions. These neutrinos are in the energy sensitivity 
range of Icecube and its DeepCore |24fE5| sub-array. 

In ^n] we brieffy introduce the astrophysical model 
and present the method of neutrino emission calculation. 
In §III| we present the results for the expected neutrino 
ffuxes and muon events, as well as the detection prospects 
with Deep Core and IceCube. A discussion and summary 
of the results is given in §IV[ 
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II. THE NUCLEAR COLLISION SCENARIOS 

We consider two types of nuclear collision scenarios, 
one where the collisions occur as a result of longitudinal 
(radial) velocity drifts, and another where they occur as 
result of transverse (relative to the jet axis) drifts of neu- 
trons from an outer jet sheath to an inner jet core where 
the bulk Lorentz factor is different. We consider both 
of these cases in the context of magnetically dominated 
jet dynamics, which differs from the usually considered 
baryonically dominated dynamics. 



A. Neutrinos From Longitudinal Nuclear Collisions 

In a magnetized outflow the bulk Lorentz factor accel- 
erates initially as 



r = 




(1) 



where Vsat = if' fa is the saturation radius beyond which 
the jet material starts to coast. (This is in contrast to 
the baryonic dominated dynamics, where F cx r up to an 
rsat ~ VTq)- 

At the radius where the co-moving baryon collision 
timescale becomes longer than the adiabatic expansion 
timescale, the protons decouple from the neutrons, and 
they continue to accelerate as F cx r^^^. The condition 
above is expressed as 



exp 



where 



(2) 



(3) 



is the total comoving baryon density in the flow, Ltot 
is the jet total isotropic equivalent luminosity, and 77 — 
Ltot / Mc^ is the dimensionless entropy or energy to mass 
outflow ratio. 

The neutrons thereafter coast with a bulk Lorentz fac- 
tor F„ — {rd/rf)^/'^ . where 



(4) 



is the decoupling radius « defined from the con- 
dition ^ (see also MRU). Here 

j/tt is a dimensionless 

parameter given by 



^ ATTcmpC^rQ 



1/6 



L33 X lO^L^f r"^/' 



(5) 



where cttt ^ 3 x lO^^^cm^^. Decoupling occurs at < 
fsat if > flTT, otherwise, if the condition ^ is met at 
a radius above rsat j decoupling never happens and the p 
and n just coast together. 



Beyond the decoupling radius, the accelerated protons 
collide longitudinally with the neutrons, which have a 
smaller Lorentz factor. The collisions becomes mostly 
inelastic when their relative Lorentz factor Trei ^ 1-3; 
see eq. ([s]). In the star frame , each shell encompassed 
within (r, r + dr) contributes to the pionization optical 
depth for those protons by an amount 



drir) — nn(r)A/3a^dr 



(6) 



where A/3 is the relative speed between protons and neu- 
trons in the lab-frame, A/3 = [1 - T-'^{r)]^/^ - [1 - 
F-2(r,)]V2 

To estimate the pion spectrum, we can, for simplicity, 
assume that for each collision, in the center of mass frame 
the proton and neutron are approximately at rest after 
the collision, with a maximum number of pions created, 
which are also approximately at rest (we will use a more 
detailed treatment in SUB I. Therefore, the invariant en- 
ergy is 



[p^ + rf)2 = V2(l + r.ei) 



(7) 



where the proton Lorentz factor viewed in neutron co- 
moving frame is 



1 Tp r„ 

^ ^ n ^ p 



(8) 



This relation is valid when both Fp , F„ ^ 1 . 

The maximum number of pions that can be created 
(the pion multiplicity) is 



A, = (V2(l + r,,,)-2)/(m,/mp) (9) 
with a Lorentz factor 

rp(r)+Fd 



r.(0 



(10) 



The radius-dependent probability for a proton to interact 
with a neutron is 



P{r) = e^^('') 



(11) 



where r(r) — J^^ dr and drir) is given by eqn. (|6| The 
number of pions created per proton is 



(12) 



where the maximum interaction radius is estimated as 
rmax ^ ctdurationi AP ~ sat (the last equaHty is for 
the choice of nominal parameters in MRU. Otherwise 

1'raax ^ ^^n—decay^n)- 



^ neglecting random thermal motions and transverse collisions, 
which are discussed in the next section. 
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Beyond the decoupling radius r^, the protons continue 
to be accelerated until reaching rgat- When the energy 
is above the threshold to create one pion in a collision 
with a neutron, we define this radius as, e.g. ri . At a 
larger radius, protons with greater energies create more 
pions per collision. The pion Lorentz factor as a function 
of their production radius is given by eqn. (10). The 



inelastic interactions starts from ri and may last to a ra- 
dius above r^at- Beyond r^at, the protons coast with a 
Lorentz factor rj and (provided rj > tItt) create pions with 
a monochromatic energy (in the first-order approxima- 
tion used in this section). In fact, the resulting neutrino 
spectrum, after the pions decay, will be broadened due 
to various factors, e.g. a) thermal motions of both pro- 
tons and neutrons; b) energy dispersion of the created 
charged pions; c) kinematics of the pion and muon de- 
cay process; d) pion and muon cooling and re-scattering. 
In this §11 A| the above four factors are not considered 
in the calculation, except in a qualitative way. This is 
adequate because, as discussed in illB and in Fig.[l], 
it turns out that the dominant process for neutrino pro- 
duction is through transverse nuclear collisions. As an 
example, for a case with rj = 500 , tq = roj , Ltot = -^54 
, Lp = L„ = O.bLh (where the sub-index p, n, b refer 
to proton, neutron, baryon luminosity respectively), the 
charged pion spectrum from longitudinal nuclear collision 
is shown in Fig.[l]. For simplicity reasons, a Gaussian 
distribution is assumed to represent the dispersion. The 
pion spectrum from transverse collisions is also shown in 
the figure, anticipating the results of the more detailed 
treatment of the latter in §IIB| and after. 



where < Trei > is the average relative Lorentz factor 
between the neutron and the baryons encountered along 
its path We estimate the number of neutrons which 
have drifted into the core in a timescale t± (all measured 
in the star frame) as 



Nj_{tj_) - vrr^ 



K{t)dt 



(14) 



where the diffusive fiux (/)n,(t) is from eqn. (14) in MRU. 
Meanwhile, the number of baryons (n and p) passing lon- 
gitudinally through the jet core is 



(15) 



The number of collisions per baryon in the core is 
roughly N±/N^^, because the pionization optical depth 
is T ^ 1 for those neutrons. (The case N±/N^i^ ^ 1 is 
less likely for typical jet parameters, and would involve 
protons in the core undergoing multiple scatterings re- 
sulting in cascades, which is beyond the scope of this 
paper). From the equations above, 

iV^/iV|| ^ 2ry-J,r;-3(i _ T'~Y^^iyrj/t^y/W (16) 

where r± is expressed in units of ro and T' is the co- 
moving temperature in units of nipC^ (T' « o(l) in this 
case, see MRU). The parameter y is the ratio of neu- 
tron density in the sheath to baryon density in the core, 
y = rin.out/'nb, in the star frame. 



B. Neutrinos From Transverse Nuclear Collisions 



In more realistic jet models, the jet properties vary 
in the transverse {9) direction. Hydrodynamical simula- 
tions indicate a smaller Lorentz factor in the jet outer 
regions (jet edge) than that in the jet core [30], and 
qualitatively similar results appear also in some MHD 
outfiows |31j. As a simple ideal model, we consider the 
transverse structure of the jet in the region outside the 
star to be represented by a two-step function, consisting 
of an inner jet core with F give by eqn. ([T|) and a slower 
outer jet, or jet sheath, with a saturation Lorentz factor 
of riout = 10^??otit.2- Both inner and outer jets will have 
been populated with protons and free neutrons already 
near the black hole, where any nuclei present would have 
been photo-dissociated. Due to thermal diffusion, neu- 
trons in the jet sheath can drift sideways into the jet core 
and interact with baryons in the core. For significant ef- 
fects, this requires the neutrons in the jet sheath to drift 
through a substantial transverse distance ^ rO into the 
core. This condition can be roughly estimated as 

> r^le < Trei > ro/r/, (13) 



C. Parameters of the Model 

The neutrino spectrum at Earth depends on a number 
of parameters. Among these are the total luminosity Ltot 
of the jet, which is initially mostly in magnetic form; the 
"baryon" luminosity Lf,, which is the the dominant energy 
form beyond the saturation radius Vgat — V'^'i^o (the ratio 
of these two being defined as £{, = Lb/Ltot)', the pho- 
ton luminosity produced by dissipation process around 
the photosphere Lph (whose ratio to the total initial lu- 
minosity is defined as e^/j = Lp^/ Ltot)- In the inner jet 
core, the proton and neutron luminosity are assumed, for 
simplicity, to be Lp — L„ — {l/2)Lh throughout the pa- 
per We take for the ratio y of baryon density in the 



^ The path calculation is a well-defined but complicated problem. 
To derive eqn. | |13[ we set the pionization optical depth = 1 
along the neutron path for which its transverse displacement 
amounts to rO. The path itself is not transverse to the jet axis, 
due to relativistic beaming according to different sheath and core 
bulk velocities. < F^e; > depends on the details of the jet trans- 
verse structure, and the density along the neutron path varies 
due to the jet dynamics. However, eqn. l |13| serves as a rough 
estimate. 

^ The outer jet neutrons can also collide with jet core neutrons 
and produce neutrinos. However, the relative Lorentz factor with 
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outer jet to that in the inner jet a nominal value y = 0.01. 

and riout 

3|4 



The discussions on results from different rjcore 



are presented in Figure 
opening angle of — O.OT 



We adopt a nominal jet 
other values being discussed 



in Figure [s] . We also adopt as a standard burst duration 
in the source frame a value of t = 20s, which is a rough 
average value for long bursts. The redshift-distance rela- 
tion used is that for a standard ACDM cosmology. 



III. GEV NEUTRINOS AND THEIR 
DETECTABILITY IN DEEP CORE AND 
ICECUBE 

A. Neutrinos from a single GRB 

The neutrinos are produced by nuclear collisions be- 
tween protons and neutrons leading to pions, the charged 
pions subsequently decaying]^ as 



(17) 



The neutrino flavor mix produced at the source is deter- 
mined by eqn. 17 , but as a result of neutrino oscillations, 



the neutrino flavor received at Earth depends on energy 
and distance. Since in this paper we discuss the general 
case, not a speciflc GRB, we approximate the received 
neutrino flux as having equal numbers in all three fla- 
vors. 

To calculate the neutrino spectrum from nuclear colli- 
sions the use of a numerical code is desirable. The com- 
monly used PHYTIA-8 code uses a minimum threshold 
energy of Ecm = 10 GeV, which for the energies consid- 
ered here leads to inaccuracies. For this reason, we have 
used two different numerical methods. One method uses 
the publicly available code of |32j. These authors use a 
parametrization formula for 7,e^,i/,and P which is carried 
out separately for diffractive, non-diffractive processes 
and resonance-excitation processes, with a logarithmic 
rising pp inelastic cross section with Tp. The secondary 
particle spectra are initially extracted out of events gen- 
erated for monoenergetic protons (0.48806^ < Tp < 
512TeV), using several simulation programs. The spec- 
tra are then fltted by a common parametrized function, 
separately for the physical processes listed above. Fi- 
nally, the parameters determined for monoenergetic pro- 
tons are fitted as a function of proton energy. The proce- 
dure was repeated for all those secondary particle types 



neutrons is smaller than that with the accelerated protons in the 
jet core, so the neutrino production through nn collisions is less 
efficient than for pn collisions. Only the latter are discussed in 
this paper. 

* The charged Kaon leading decay channel is similar and can be 
approximated as effective pions. However, the number ratio of 
Kaons to pions produced by hadron collision at these energies is 
less than 0.1 and can be neglected for the purpose of this paper. 



mentioned above. In order to approximate better the ex- 
perimental data at lower energies, two baryon resonance 
contributions have been included, one representing the 
A(1232) and the other representing multiple resonances 
around 1600 MeV/c^'^. However, as pointed out also by 
|32| , the pion mean energy in this code is slightly under- 
estimated at incident proton kinetic energy (in the fixed 
target lab frame) of Tp ^ 2 GeV and above, compared 
to experimental data. Therefore in this paper we have 
corrected this discrepancy by multiplying the resultant 
neutrino energy by a factor of 1.3 using an approximate 
fit to the Fig. 5 of their paper. This causes a slightly 
over-estimated pion mean energy near the threshold en- 
ergy Tp ^ 0.3 GeV; however, this energy range is not of 
interest for our purposes in this paper, since it results in 
a non-detectable neutrino fiux is associated at the asso- 
ciated energy. We refer to this numerical calculation as 
method A. 

We have also developed a different code, which is inde- 
pendent of method A, in order to cross-check the validity 
of method A around energies Tp > 4 GeV, and in order 
to gain better transparency on the underlying physical 
processes where this is not otherwise made explicit in 
method A. We refer to this second method as method B. 
The energy Tp > 4 GeV is of particular importance here 
for at least three reasons: i) it arises naturally in the as- 
trophysical model considered here; ii) nuclei in this range 
can produce substantial neutrino fluxes; iii) Icecube and 
its DeepCore sub-array neutrino detectors are sensitive to 
the details of the neutrino spectrum in the energy range 
10—1000 GeV. In method B, the charged pion spectrum is 
approximated by a radial scaling based on the appar- 
ent Feynman scaling violation at xu = E/{y/s/2) <C 1. 
The meson decay kinematics are well established, and in 
method B we follow the formulation of [34 

A comparison of the muon neutrino and anti-neutrino 
spectra at the source calculated using method A and B is 
plotted in Fig.|2|, for an incident proton energy Tp = 3.8 
GeV. It is seen that the two results agree well with each 
other. We note that a) both methods would result in a 
very small number of of neutrinos which violate energy- 
momentum conservation. Although the total amount of 
energy involved in these neutrinos is well below a fraction 
^ 10"'^, and thus negligible, we have nonetheless applied 
a cutoff in the spectrum beyond the energy where this 
occurs for both method A and B. b) both these methods 
simulate pp collisions, instead oipn collisions. For pp col- 
lisions, the TT"*" multiplicity is therefore greater than that 
of TT~ due to charge conservation near the pion creation 
threshold energy. At y/s ^ few GeV or higher, the two 
multiplicities tend to equal each other. In method A we 
sum the neutrino and anti-neutrino from both 7r+ and tt^ 
channels. In method B, Hillas used one fitted formula to 
describe both tt"*" and tt" spectra. Thus, for the summed 



^ Noting a typo in this reference, where in his Eq.[14], second line, 
rji, should be replaced by ^. 
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at low Lp, Ln grows fast with it and then sta- 
10^^ erg/s. This effect is mainly due 



neutrino and anti-neutrino spectrum resulting from the 
decay products of these mesons, the discrepancy between 
pp and pn collisions becomes less noticeable. 

The muon neutrino and anti-neutrino differential num- 
ber fluence at Earth for a single burst (neutrinos per en- 
ergy decade, integrated over the duration of the outburst) 
is shown, after oscillations, in Figs. |3|4|5 . In Figure 
3p , we note a saturation effect in the dependence of L^, 
on L 

bilizes above L 

to the change of relative Lorentz factors between 77, 
and rjout by choices of different astrophysical parameters. 
Trei grows with Lp if we fix other parameters, and satu- 
rates at some Lp^o level (may be already saturated at the 
low end of Lp in parts of our parameter space. A higher 
Trei is associated with more secondary leptons (including 
neutrinos) per pn collision. However when saturated, the 
received neutrino flux grows at a slower rate with Lp, be- 
cause the increase of fluence is then only due to the fact 
that we have more protons in the jet. 

As discussed in §2, most of these neutrinos come from 
transverse drift collisions, rather than from radial drifts. 
We have used parameters representative of standard long 
GRB, extending from moderate to high intrinsic lumi- 
nosity, and the fluence is calculated for a nominal red- 
shift of z = 0.1, corresponding to a luminosity distance 
of 450 Mpc in a standard fi„j = 0.28, = 0.72, 
h = 0.72 cosmology. This is at the lower end of the 
classical redshift distribution, since for the 5-100 GeV 
neutrino energies considered here only rare, very nearby 
bursts might be expected to be detectable individually, 
due to the detector effective area decrease with energy. 
The typical neutrino fluences in Figs. 3|4|5 
2 ~ 1.6 X 10 " " " " ' ' 



are < GeV 



cm"" ~ i.t) X iU~'^ erg cm~ for z = 0.1. Since pp col- 
lisions are expected to produce a comparable amount of 
energy in photons from tt" decay as in neutrinos, it is im- 
portant to check that such a photon luminosityjjdoes not 
violate electromagnetic observation constraints. There 
are no gamma-ray detections at TeV energies so far, and 
the smattering of GeV detections involve bursts typically 
at much higher redshifts than considered here, so the pos- 
sible constraints are mainly the 20-300 keV energy range 
fluence statistics from the BATSE 4B |33 compilation 
of bursts. We can conservatively assume that at most a 
fraction < 0.3 of the total 7r° photon energy will appear 
in the 20-300 keV range, i.e. < 5 x 10~* erg cm~^. Such 
fluences are marginally compatible with the BATSE 4B 
statistics, which does include some objects with fluences 
< 10"^ erg/cm^ (BATSE does not provide redshift infor- 
mation, but it is known that most BATSE bursts are at 
z > 1; here we considered our bursts at z < 0.1, and if 
these were placed at the typical z > 1, they would show 
a BATSE fluence < 5 x 10~^ erg cm~^). Thus, our elec- 



tromagnetic fluences appear compatible with the current 
observations. 

The number of muon events arising from the above 
incident muon neutrino and anti-neutrino fluences is cal- 
culated based on the speciflc detector characteristics. For 
the Icecube full 86-string operation, the effective area is 
larger than that of Deep Core at ~ 100 GeV, these ef- 
fective areas being given by j^S]. For a detection rely- 
ing exclusively on neutrinos one would need to consider 
GRBs which show at least one muon event. However, 
bursts able to give an average of > 0.01 muon events 
are also of interest, because of the expectation of natu- 
ral fluctuations in the distance or in the luminosity, and 
because in some cases temporally coincident electromag- 
netic observations can be expected. The burst values of 
L-y and redshift z yielding different muon event numbers 
for different burst parameters are shown by the contours 
in Fig.|6]. 

The average rate per year at which GRBs occur pro- 
ducing > 1 muon events can be estimated using the 
GRB luminosity and redshift distributions discussed in 
pil B[ For a baryon to photon luminosity ratio of 10, i.e. 
eb = lOee ~ 1, with inner to outer jet Lorentz factor 
contrasts of ricore = 10??o«t = 300, 700, 1000, respectively, 
this average rate is expected to be around 0.0006/yr, 
0.04/yr and 0.06/yr. 



B. Diffuse Neutrino Flux 

The diffuse neutrino flux from all GRBs in the sky can 
be calculated using a GRB luminosity distribution (lumi- 
nosity function, LF) and a redshift distribution (RD)Q 
Here we adopt for these the functions given in |37| , 



(f>{Ly) CX 



{L^/L,Y'^ Lr, 



Rgrb{z) 



{i + zy 
[i + zY 



Z < Zl 
Z > Zl 



(18) 



(19) 



Here equation |18| is the luminosity function, and Eqn. 



19 is the redshift distribution function, L^ is the peak 
pEoton luminosity (typically in the 0.1-0.3 MeV range), 

Lrmn = 10^° Crg/s, L, = lO'^^.S g^g/g^ ^^^^ ^ ^q54 

erg/s, nil = —0.2 , m2 — —1.4 , rii = 1.0 , n2 = —1.4 
, Zl = 3 , and we have used these values from the |37| 
parameter ranges which best reproduce the actual GRBs 
with measured L^ and z statistics, e.g. Fig. 2 of |37| . 
Fig. 4 of |i3S], or [35] • We note that, especially in the 
range z < 0.5, the index of the redshift distribution, i.e. 
the rate, is very uncertain, due to poor statistics. The 



^ An evaluation of the final photon spectrum would require the use 
of a detailed electromagnetic cascade code, which is beyond the 
scope of this paper. 



^ For the method of calculation, see e.g. Appendix B of I36| 
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differential co-moving rate of GRBs at a redshift z is 
given by 



R{z) 



Rgrb{z) dV[z) 
(1 + z) dz 



(20) 



where V{z) is the comoving volume in the ACDM cos- 
mology model adopted, with H.^ = 0.28, fiy = 0.72 and 
Ho = 72 Km/s/Mpc. The differential number of GRB 
per unit redshifts is given by 



dN = pQ(f){L^)R{z)dlogLjdz 



(21) 



where po is the normalization factor. In this paper, we 
normalize the total electromagnetically detecte d GRB 
rate to 300/yr in the range given by Eqns. 18|19 . 
The neutrino flux depends on the baryon luminosity Lb, 
and the ratio between Li, and adopted in JIII A is 
Lh/Ly « 10. This ratio is characteristic of many hadronic 
GRB models, the implied photon radiative efficiency of 
10% being moderate. The resulting diffuse neutrino 
fluxes (neutrinos per year) are shown in Fig.|7], and are 
discussed in 3IV| 



IV. DISCUSSION 

We have calculated the neutrino emission in the range 
of a few GeV to a few hundred GeV arising in magnetized 
coUisional GRB models, such as have been recently used 
for interpreting the electromagnetic properties of these 
objects. The neutrino emission considered here arises 
partly from longitudinal proton and neutron collisions 
following 'their decoupling in the jet, and in larger part 
from collisions caused by transverse thermal drift of neu- 
trons from an outer jet sheath into the jet core, hav- 
ing different bulk Lorentz factors. This neutrino produc- 
tion model differs from the commonly considered photo- 
hadronic pj models, and it also differs from previous pn 
model calculations in incorporating explicitly the (domi- 
nant) transverse drift collision effects. Also, the emission 
region characteristics are here determined by the magnet- 
ically dominated dynamics, differing from those in pre- 
vious neutrino calculations, which mostly use baryonic 
dynamics. Furthermore, the pn neutrino spectra are cal- 
culated numerically, using two different codes which are 
suited for the GeV range energies considered. 

Our present results indicate that for the burst parame- 
ters suitable to explain the photon spectra and the MeV- 
GeV photon lags indicated by the recent Fermi satel- 
lite observations, a low level of neutrino emission is ex- 
pected at ~ 5 — 100 GeV energies. This is in the sensi- 
tivity range of the Deep Core sub-array of IceCube, and 
extends into the lower range of the main IceCube array. 

For a neutrino detection of an individual burst, un- 
aided by a coincident electromagnetic detection, one 
would require > 1 muon event. For neutrino-induced 
muon track events in IceCube and Deep Core, consider- 
ing the angular resolution (e.g. D, < 10^ deg^) and time 



search bins ~ o(l) x 20s per burst, the atmospheric neu- 
trino background is low. During the total sum of the 
search bins (e.g. assuming 300 electromagnetically ob- 
served long GRBs in a year, from the atmospheric neu- 
trino spectrum the expected number of background muon 
events is < o(l) x 0.09. 

However, any GRBs which might be expected to yield 
> 1 muon event need to be at the high end of the 
luminosity function, and located at very low redshifts, 
L-y/iso ^ 10^'^ erg/s and z < 0.1. From the discussion 
of §111 B| the occurrence of such GRB is estimated to be 
very rare, < 1/17 per year. This estimate is uncertain, 
because of the poor statistics in the determination of the 
redshift distribution in this range. Such bursts would 
in general also be detectable by photon detectors such as 
Swift and Fermi, except for Earth occultations or possible 
outages. For the more frequent weaker or more distant 
bursts, taking into account fluctuations in the average 
quantities, a neutrino observation correlated with a pho- 
ton detection can narrow the time bin search, increasing 
the effective sensitivity of the detection. The rate of oc- 
currence of such lower fluence GRBs can be calculated 
from the luminosity function and redshift distribution, 
and is shown in Fig. [6]). 

Because of the low occurrence rate of bursts which can 
be expected to be detected individually, it is useful to 
consider also the diffuse neutrino fluxes. The same lumi- 
nosity function and redshift distribution as above were 
used for this, as discussed in §111 B| extending the inte- 
gration to all bursts with z > 0.01. In Fig.[7] we show 
the diffuse neutrino fluences over the period of a year, for 
a set of burst parameters rjcore = 10?7out = 300, 700, 1000 
and £(, = lOee ~ 1. For these paramaters, the correspond- 
ing number of average expected muon events in Icecube 
and its DeepCore array are estimated to be 0.03, 0.4 and 
0.5 per year over the whole sky, respectively. Due to the 
poor statistics in the rate of low redshift GRBs, the un- 
certainty in these numbers could be a factor ~ 2. The 
Icecube IC40+IC59 muon event rates per year being used 
to set upper limits on the TeV-PeV neutrino flux from pu- 
tative Waxman-Bahcall GRB models (different from the 
present model) are larger than the the event rates dis- 
cussed here, but within one order of magnitude of those 
for the highest 77 values. Thus, imposing weak limits may 
perhaps be possible in the long term. 

In conclusion, both the individual burst fluences and 
the expected diffuse flux in the 5-100 GeV range are sig- 
nificantly low. Fortuitous fiuctuations above the mean 
values could increase this somewhat, but any conclu- 
sions based on the current IceCube and Deep Core ar- 
rays are likely to require years of data accumulation with 
the full array. The proposed upgrades to these installa- 
tions would help, but a next generation of larger effective 
volume neutrino detectors could be the best way to ac- 
celerate the detection or non-detection of GRB neutrinos 
in this energy range, and to test GRB models such as 
discussed here. 
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100 150 200 250 300 



FIG. 1. Charged pion spectrum at the source from longitu- 
dinal and transverse nuclear collisions (analytical approxima- 
tion), normalized to the Earth observer frame. The charged 
pions decay in the source, and are observed at Earth only via 
their neutrino decay products. The component LI is from 
longitudinal pn collisions at radii rd < r < Vsat, while the 
component L2 arises at r > Vsat- The component T is from 
transverse drift pn collisions (see 511 B I. A Gaussian disper- 



sion is assumed for simplicity to represent the broadening 
of the real spectrum arising from various effects (see ill A I. 
The dominant pion (and neutrino) production comes from 
the transverse pn collisions, as discussed in §11 B| We assume 
a source with = O.lLtot = lO^^erg/s, rj = 500, Ojet = 0.01 
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FIG. 2. The sum of the muon neutrino and anti-neutrino spec- 
trum from nuclear collisions, calculated in the stellar frame 
using the numerical method A (curve marked by " + ") and 
the method B (curve marked by "o"), as discussed in |1II A 
Both spectra are normalized to one pp collision event. In the 
outer jet comoving frame (where target protons are at rest) 
the incident proton kinetic energy is Tp — 3.8 GeV. In the 
case shown here the outer jet Lorentz factor is rjout = 70 and 
the inner jet's is r^core = 700. 



and z = 0.1 (see iUC\ 
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FIG. 3. The u^+i?^ fluence in the Earth observer frame from a 
single burst, integrated over the outburst duration, for a nom- 
inal redshift oi z = 0.1. The other nominal parameters are the 
same as in Fig.[l], except for rjcore ~ Wrjout ~ 300, 700, 1000, 
indicated with dashed, solid and dot-dashed lines respectively. 
The lines in each style are arranged from top to bottom in the 
sequence Ltot ~ Li, = 10^^ 10^^ ■^ 10^^ 10^^=* ■^ 10^^ erg/s. 




FIG. 4. The //^ + fluence in the Earth observer 
frame from z — 0.1, for different parameters. Top 
panel: solid lines from bottom to top are for rjcore = 
300, 400, 500, 600, 700; dashed lines from bottom to top are 
for ??core = 2000, 1500, 1000, 900, 800, with a fixed rjout = 100 
Bottom panel: solid lines, bottom to top are rjcore = 
300, 400, 500, 600, 700; dashed hues, bottom to top are rjcore = 
1500,1200,1000,900,800 with a running rjont = O.lrjcore- 
Both panels are for Lp — 10^* erg/s, the value of ijcore yield- 
ing the largest neutrino flux being around rj ~ 700. When 
rjcore approaches 300 the pion multiplicity becomes smaller 
due to a smaller relative Lorentz factor between outer and 
inner jet particles (top figu re. Eg. [9]) or due to a slightly 
smaller diffusion rate (Eq. |l6| ). For very high rjcore ^ 10^, 
the pion multiplicity similarly decreases. 
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FIG. 5. The Vfj, + number fluence in the Earth observer 
frame from z — 0.1, calculated for rjcore ~ Wrjout = 700, 
Lt ^ Lb = for various jet opening angles 9 jet = 

0.01, 0.02, 0.03, 0.04, 0.05, 0.1 (top to bottom). 
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FIG. 6. Expected number of muon events in the full 86- 
strings IceCube and its Deep Core sub-array, from a sin- 
gle GRB with Tjcore = Wrjout = 300 (top panel), 700 (middle 
panel), 1000(bottom panel), from various redshifts z and for 
different photon luminosities L^, assuming a baryon to pho- 
ton luminosity ratio of 10, or ei, = 2ep ~ lOe-y. The effective 
detection areas are taken from [25!, with the angular position 
averaged over the northern sky (the effective areas have some 
dependence on the incident angle of the neutrinos). These 
contours show the and z ranges that give 0.01 — 2 muon 
events 
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FIG. 7. The diffuse muon neutrino and anti-neutrino fluence 
per year (after oscillations) for all GRBs down to z > 0.01, as 
discussed in §IV| usin g the luminosity function and redshift 
distribution of §111 B] The conventions and parameters are 
the same as in Fig.|6], except for using ?7core ~ Wriout ~ 300 
(dashed), 700 (solid) and 1000 (dot-dashed). 



